The effect of random coupling coefficients on decoherence 
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The aim of this letter is to analyze the effect on decoherence of the randomness of the coupling 
coefficients involved in the interaction Hamiltonian. By studying the spin-bath model with computer 
simulations, we show that such randomness greatly improves the "efficiency" of decoherence and, 
then, its physical meaning deserves to be considered. 
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Introduction. Environment induced decoherence 
(EID) is usually tested by means of computer simula- 
tions on models whose coefficients are, in general, taken 
as random with the idea of simulating a generic situation 
(see, as examples, [l|, 0, 0: 0], @, @)- This widely 
spread strategy is reasonable for the coefficients of the 
particles' states, but it is not natural for the coupling co- 
efficients involved in the interaction Hamiltonian: there 
is no reason to suppose that particles of the same nature 
interact with each other with different strengths. The 
point would be irrelevant if the randomness of the cou- 
pling coefficients did not affect the overall phenomenon. 
The aim of this letter is to show that this is not the case: 
the randomness of the coupling coefficients greatly im- 
proves the "efficiency" of decoherence. This result has to 
be taken as a warning against the uncritical use of ran- 
dom coupling coefficients for drawing conclusions about 
decoherence. We will argue for this claim by means of 
the analysis of a well-known model. 

The spin-bath model. This is a very simple model 
that has been exactly solved in previous papers (see [l|). 
We will study it from the general theoretical framework 
for decoherence presented in a previous work [7] . Let us 
consider a closed system U = S + E where (i) the system 
S is a spin- 1/2 particle P represented in the Hilbert space 
Hs, and (ii) the environment E is composed of N spin- 
1/2 particles Pi, each one represented in its own Hilbert 
space Hi- The complete Hilbert space of the composite 

/ N \ 

system U is 7i = TLs ® I ® I ■ in the particle P, 



the two eigenstates of the spin operator S s V in direc- 
tion V are Iff) and |JJ.), such that S s ^ Iff) = \ |fr) and 
&s if = — \ W)- ^ n eacn particle Pi, the two eigen- 
states of the corresponding spin operator S i in direc- 
tion V are |Ti) and such that S t |T») = \ \\i) and 
l-U) = — | Therefore, a pure initial state of U 
reads 



IV>o> = IfT) + b W» ® 6d(ai\ U) + Pi\ U)) (1) 



where the coefficients a, b, cti, /3 i are such that satisfy 
a| 2 + |6| 2 = 1 and ai| 2 + |/3j| 2 = 1. The sclf-Hamiltonians 
Hs and He of S and E, respectively, are taken to be zero. 



Then the total Hamiltonian H = Hs + He + Hse of the 
composite system U results (see [l[, 0) 
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H = H SE = S s ^ ® 2 9iSi^ ® (g) h 
i=l \ j^i 



(2) 



where Ij is the identity operator on the sub- 
space Hj, S s ^ = ± (|fr) (fr| - |^> <^|) and = 
3 (Itt) (Til - \U) CUD- Under the action of H = H SE , the 
state \yb Q ) evolves as \ip(t)) = a |fr) |£ t (f)) +b\i^) |£ 4 (t)) 
where \£-ff(t)) = |£|l(— i)) and 



\£ t (t)) 



N 

®(a z e-^ 2 \U)+(3 ie ^ 2 II,}) (3) 
i=i 



An observable of U, O G = H <E)H, reads 



O = 



+s U I1r> 

+«w I4> (frl 
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+e | ||.) 
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(4) 



In the typical situation studied by the EID approach, the 
system of interest S is simply the particle P. Therefore, 
the relevant observables Or are obtained from eq.(j4|) by 
making ei» = ef) = 1 and e{ 1 ) — 0: 
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(>n=\ wW*'! ® =Os®J B (5) 

^s,s'=ir,tt / \i=l / 

The expectation value of these observables in the state 
is given by 

(Ofl)^(t) = M 2 S M + l^l 2 S M + 2 Re[a6* s M r(f)] (6) 
where (see 0) 

r(t) = = JJ [la^e-^* + l/J^e^ 4 ] (7) 

i=l 

It is clear that the time-behavior of {Or)^) depends on 
the behavior of r{t). In order to know such a behavior, 



we performed numerical simulations for 

N 



<t)\ 2 = n(i^i 4 +i^i 4 + 2 i^i 2 i^i 2cos2 ^) 
1=1 

N 

= UMt) (8) 



where the la; 



generator, and the were computed as I/3J 2 



were obtained from a random-number 

1 - 

\ai\ z . But before presenting the computer simulations, 
we will study the Poincare time for this model. 

Analyzing the Poincare time. Each /»(£) of eq.® 
comes back to its initial value for the first time at a time 
tpi, such that 2gitpi — 2ir ==> tpi — n/gi. In turn, |r(i)| 2 
comes back to its initial value when all the fiit) do it. 
Therefore, the Poincare time tp of this model is the time 
when all the fi{t) come back to their initial values for the 
first time. Let us consider three cases: 

(a) All the gt have the same value: gi = g, for all i. 
So, the mean value is gi = g. In this case, all the 
fiit) come back to their initial values at the same 
time tpi — n/g. Therefore, the Poincare time is 
tp = n/g = Tr/g~i, and it does not depend on the 
number of particles N. 

(b) All the gi are such that g^ = nig m { n , with Hi £ N. 
In this case, tp is the largest tpi, corresponding to 
the smallest g i} g min : t P = n/g mia . So, given a 
<?minj the Poincare time tp does not depend on N. 
Since the mean value is ~gl > g mm , then tp > ir/gi. 

(c) All the gi are random. Since in the computations 
the gi are rational numbers, we can express them 
as gi = Pi/ qi, with ft, % £ N. If we make tp = -kQ, 
the number Q has to be such that Q = riiqi/pi for 
all i, with m G N. Then, m = Qpi/qi- Since the 
Pi and qi are random natural numbers, the least Q 
that guarantees that ni is a natural number for all 
i is Q = rL=i 1i- Therefore, 
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t P = ttQ = tt]^[' 



(9) 



In turn, Q is larger than any qi, more larger as N 
increases. Then, for N large and for any g+ = Pi/qi, 
the Poincare time is tp = irQ 3> Ttqi > nqi/pi — 
ir/gi = tpi, and also tp ir/gl. Moreover, 
the order to magnitude of can be estimated as 
tp = ttQ ~ ^qJ N (^/9i) N ' when the coupling 
coefficients are random, the Poincare time increases 
exponentially with the number of particles. 

Case 1: Homogeneous environment. Let us 

consider the case where all the particles of the environ- 
ment are of the same kind. In this case, the reasonable 
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FIG. 1: Plot of \r{t)\ 2 given by eq.fl8|, for N = 100 and 
9i = 9 = 0-5. 
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FIG. 2: Plot of \r{t)\ 2 given by eq.flSJ, for N = 100 and 
random gi g [0.4,0.6]. 



assumption is that the particle P interacts in the same 
way with all the environmental particles and, as a con- 
sequence, gi — g for all i (Case (a)). The time-behavior 
of \r(t)\ 2 of eq.©, for N = 100 and g. t = g = 0.5, is 
plotted in Figure 1: we can see that the Poincare time 
tp = n/g = 7r/0.5 ~ 6.28 is not sufficiently larger than 
the decoherence time to consider the result as an effective 
decoherence that may lead to classicality. 

Now let us consider the case that P does not interact 
in the same way with all the environmental particles. In 
particular, we will suppose that the coupling coefficients 
are random in an interval \gl — Ag,g~i + Ag] around the 
mean value gi. The time-behavior of |r(£)| 2 of eq.©, for 
N = 100 and gi random, with gi — 0.5 and Ag = 0.1, 
is plotted in Figure 2. In this case, the Poincare time is 
much larger than {-n /~gi) N = (6.28) 100 (Case (c)), a value 
that can be considered infinite for all practical purposes. 
Therefore, as Figure 2 shows, it can be legitimately said 
that the particle P decoheres and may become classical. 

Case 2: Non-homogeneous environment. In 
this case we will consider an environment whose particles 
are not all of the same kind, and the particles of each kind 
j interact with P through their own coupling coefficient 
gy. Ni particles through g\, N2 particles through 92, 



and N p particles through g p , such that 



p 



N 4 



N. 



Then, the |r(£)| 2 of eq.© can be rewritten as: 
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FIG. 3: Plot of |r(t)| 2 given by eq.([ID]l. for N = 100, Ni = 91, 
^2=^3=^4 = 3, gi = 2.4, ff2 = 1.2, g 3 = 0.6 and g 4 = 0.3. 



FIG. 4: Plot of \r(t)\ 2 given by eq.flTT}, for TV = 100, Ni = 
91, ^2=^3=^4 = 3, and random g u G [1.6,3.2], g 2i G 
[0.8,1.6], g 3l G [0.4,0.8], and g u G [0.2,0.4]. 
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\r(t)\ 2 = 1] /<(*) = ( 10 ) 
i=l 

/ Ni \ / AT 2 \ / 



i=l 



>*=1 



where each particular product is the contribution of each 
kind of particles. In particular, we will consider a situ- 
ation where the environment is composed of N particles 
that are almost all of the same kind, with the exception 
of a slight "contamination" of particles of different kinds. 
The time-behavior of \r{t)\ 2 of eq.lJTDJ), for N = 100, 
JVx = 91, N 2 = N 3 = N 4 = 3, .gi = 2.4, g 2 = 1.2, g 3 = 0.6 
and (?4 = 0.3, is plotted in Figure 3. Since all the gi are 
such that gi = nig min — rii g 4 , this situation corresponds 
to Case (b). Then, the Poincare time can be easily com- 
puted as tp — 7r/g m j n = 7r/0.3 ~ 10.43. Therefore, in 
Figure 3, the peak in 10.43 is the Poincare time, but the 
peaks around it are not due to the recurrence of \r(t)\ 2 . 
This means that the particle P in interaction with this 
"contaminated" environment does not decohere. 

Now we will consider that P does not interact in 
the same way with the particles of a same kind j, but 
the coupling coefficients gji are random in the intervals 
[gji — A<7j, ~g]i + Agj] around the corresponding mean 
value gj[: Ni particles with g u G [#i7 - Apr, TfiJ + Agi], 
N 2 particles with g 2l € [gii - Ag 2 , gYi + Ag 2 ] , ... and 
N p particles with g p i G \g^l — Ag p ,~g^i + Ag p ], such that 
p 

J2 N J = N - Then ' the l r (*)l 2 of e( i- : 

i=i 
as: 



can be rewritten 



N 



<t)\ 2 = n /*w = 



(ii) 



i=l 



' N 2 



\i=l / Vi=l / \i=l t 

In particular, we will study a situation similar than the 
previous one with respect to the values of N and of the 
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FIG. 5: Same as Fig. 4 with a different time scale. 



Nj, and where the gj used there become here the mean 
values ~gji. Moreover, in all the cases the Agj were se- 
lected as approximately the 30% of the corresponding 
mean value ~gjl. The time-behavior of |r(f)| 2 of eq. (fTTI) . 
for N = 100, TVi = 91, N 2 = N 3 = N A = 3, gTi = 2.4, 
Wi = 1.2, W = 0.6, gil = 0.3, A gi = 0.8, Ag 2 = 0.4, 
Ag 3 — 0.2 and Ag 4 = 0.1, is plotted in Figure 4, where 
we can see the drastic decoherence of the model. The 
same result is plotted in Figure 5 with a different time 
scale, in order to show the decoherence time. 

Conclusions. By means of computer simulations we 
have proved the great influence that the randomness of 
the coupling coefficients exerts on decoherence. Such an 
influence appears under two forms: (i) the huge increas- 
ing of the Poincare time, and (ii) the strong damping off 
of the peaks that preclude decoherence. Then, in both 
cases the change from constant coefficients to random 
coefficients greatly improves the "efficiency" of decoher- 
ence. It is worth stressing that, whereas the first effect 
may be expected, the second effect is not foreseeable: it 
is not easy to explain a priori why the randomness of the 
coupling coefficients washes off the high peaks present in 
the non-random case. 

On the other hand, a realistic environment is usually 
composed by a limited number of different kinds of par- 
ticles, and it is reasonable to suppose that the particle 
of interest P interacts with all the particles of the same 
kind with the same strength. So, the models with random 
coupling coefficients are not realistic. As we have pointed 
out in the Introduction, this would be not relevant if the 
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randomness of the coefficients had no significant effect. 
But now we know that such randomness has a dramatic 
influence on decohcrcncc, whether it is expectable or not. 
Therefore, the criterion for the selection of the coupling 
coefficients has to be carefully analyzed in each model, in 
order to avoid conclusions drawn from unphysical results. 
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